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Abstract 

Dynamics in delayed differential equations (DDEs) 
is a well studied problem mainly because DDEs arise 
in models in many areas of science including biol¬ 
ogy, physiology, population dynamics and engineer¬ 
ing. The change of nature in the solutions in the 
parameter space for a network of Phase-Locked Loop 
oscillators was studied in Symmetric bifurcation anal¬ 
ysis of synchronous states of time-delayed coupled 
Phase-Loeked Loop oseillators. Communieations in 
Nonlinear Science and Numerieal Simulation, Else¬ 
vier BV, Volume 22, Issues 13, May 2015, Pages 793- 
820, where the existence of Hopf bifurcations for both 
cases, symmetry-preserving and symmetry-breaking 
synchronization was well stablished. In this work we 
continue the analysis exploring the stability of small- 
amplitude periodic solutions emerging near Hopf bi¬ 
furcations in the Pixed-point subspace, based on the 
reduction of the infinite-dimensional space onto a 


two-dimensional center manifold. Numerical simu¬ 
lations are presented in order to confirm our anali- 
tycal results. Although we explore network dynamics 
of second-order oscillators, results are extendable to 
higher order nodes. 

INTRODUCTION 

Networks of oscillators have been studied for decades 
because their models can represent dynamics in a 
very wide range of fields, as astronomy, biology, neu¬ 
rology, economics, and the stock market, [ini ITTl 
HH 11 HI- Much of the research has focused to 
understand the influence that changes in the pa¬ 
rameter space have over the global dynamics. In 
that sense, the stability of the synchronization in 
the network have been explored from many ap¬ 
proaches [161 Ian uni ig. In this contibution, we 
study the stability of small-amplitude periodic so- 
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lutions which emerge near Hopf bifurcations of the 
symmetry-preserving solutions for a N-node PLL net¬ 
work, using the centre manifold theorem extended to 
functional differental equations, and the normal form 
in the center space. This work is mainly based on pre¬ 
vious results obtained on the existence of symmetry¬ 
preserving and symmetry-breaking Hopf bifurcations 
in a N-node second-order PLL network, by Ferruzzo 
et al. [5]. 

It is important to note that when the lag in the 
communications between nodes is considered, the or¬ 
dinary differential equation (ode) system which de¬ 
scribe the network dynamics, becomes a delayed dif¬ 
ferential equation (dde), whose solution lies, in gen¬ 
eral, in the function space, and its characteristic 
equation has infinitely many roots. This particu¬ 
lar kind of functional differential equations appear 
in many engineering problems nnmiH]. 

We consider the Full Phase model introduced in [5] 
to analyse stability of small-amplitude periodic or¬ 
bits near Hopf bifurcations emerging in the parame¬ 
ter space (/X, r) at the Fixed-point subspace, for the 
non degenerative case {K >1). It has been shown 
that these bifurcations can occur when a pair of com¬ 
plex conjugate eigenvalues crosses the imaginary axis 
in iether direction, from the left to the right and 
from the right to the left. The main approach used 
for the analysis is the decomposition of the infinite- 
dimentional space into a 2-dimensional center space 
spanned by the eigenvectors correponding to the sim¬ 
ple imaginary eigenvalues A = iiw, w > 0, and into 
an infinite-dimentional space “orthogonal” to the first 
one (the orthogonality condition will be defined be¬ 
low). We will follow closely the theory and proce¬ 
dures presented in [m Ennuis]. 


The Full-phase model 


node output phase 4>i(t), is: 


^ ^ (/>j) = 0, f = 1,..., N, 

i=i 

( 1 ) 

where: 


- r) - +sm{(j)j{t - t) + (j)i{t)). 

( 2 ) 

f : M. X M. ^ R, fi, K,t G R.+ , and N gN — {!}. 

The equilibria (j)^, in equation m , are: 

(/)+(n) = i ^arcsin + 2n7r^ 

, (3) 

4>~ (n) ~ 2 (^ ~ ^-rcsin + 2n7r^ 

n G K > \. For our analysis, we consider three 
main assumptions: 

(a) The critical eigenvalue A of the linearization 
of ([1]) at equilibria crosses the imaginary axis 
with non vanishing velocity, i.e. Re(A'(0^)) yf 0. 


(b) The purely imaginary eigenvalue A = iw is sim¬ 
ple. 

(c) The linearization of m at equilibria, has no 
eigenvalues of the form ikuj, k G Z — {1,-1}. 

The Taylor expansion of 0 at equilibria is: 


N OO 


E EI ^ ^+ 


r- 

j=i 




(4) 




where (jjjT '■= (f>j{t — T). Truncate the series up to the 
third-order term: 

E 1 


3 = 1 


i4>jT + sm2(j)^ - i {(pjr - 4>if + {4>jT + (j^if cos2(j)^ 

(5) 


The general model for a iV-node, fully-connected, 
second-order oscillator network, in terms of the i-th 
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i = here, for the sake of notation we 

changed 5/pi —>■ (pi. 

The vector field form x = G{xt,x-,r]), G : x 

X R^ — >■ R^'^, can be obtained by choosing = 
and x^'^ = (j)i, then, the restriction G\i, or, = 
G^^\xt, ; r]) gives: 


_ Ji) 


= XI "I + cos 2 (/)=^)a;i’ 


N 


N-1 


3 = 1 




(1 + cos2(j)^)xi^ — - (^Xi^ + sin2(/)- 

(^Xi^ — x^^^^ + ^ cos 2 (/)=*= 


*=!,...,iV. 


( 6 ) 


F = (/(!),...,/(^)J^, /(d = (/«,/«), /« = 0, 

X { “ ^ - 

37^3 ^ 

3=1 



In order to build the decomposition of the infinite¬ 
dimensional space, we need the adjoint operator as¬ 
sociated to the linear part of the linearization and a 
inner product, via a bilinear form. Associated to the 
linear part of ©, the formal adjoint equation is 


= Aliv)yit,v) + A'^iv)yit + Gv), (lo) 

The strongly continuous semigroup {T*{t)'ip){9) = 
{yt{'fp)){9) = y{t + d), defines the infinitesimal gen¬ 
erator: 


Following m m. we can represent the dynamics 
in ([ 6 ]) by the abstract differential equation: 

^Xticj)) = A{r])xt{(l)) + J^{xt{(l)),r]). (7) 

We define X := C([—r, 0],R^'^), the Banach space of 
continuous functions from [—r, 0 ] into R^^, equipped 
with the usual norm 

\m= sup |i?( 6 ')|, e C([-r, 0 ]), 

-r<e<0 


{A*{y)ip) 


^{0) ,O<0<T 

AoivVi’iO) + A^{r])'^'ip{T) ,0 = 0 

( 11 ) 


such that 4iT*{t)ip = A*T*{t)ij, ijj & X* := 
C([0,t],R^^). The natural inner product has the 
form ^n\ : 

M 

{x,y) =x'^{0)y{0) + J x"^ {s + T)Ar{y)y{s)ds, 

X € X and y € X*-, thus, we have [7]: 


Xt, in ([7]), lies in X and satisfies (T{t)(j))(0) = 
{xt{(j))){0) = x{t + 0), where T{t) is a semigroup of 
family of operators, 0 € [—t, 0 ], and 77 is a vector of 
parameters. The linear operator A(? 7 ) G Mat(2Af) is: 

7 47 m(^) ,-T< 0 <o 

^ I Ao(y)m+Ar(vM-T) ,0 = 0 

( 8 ) 

where ^ 0 ( 7 ?) := A^M := |g:|^± and. 


1 . A is an eigenvalue of ^( 77 ) if and only if A is and 
eigenvalue of A* ( 77 ). 

2. If (pi,..., (/Id is a basis for the eigenspace of ^( 77 ) 
and ipi,... ,ijjd is a basis for the eigenspace of 
A*{r]), construct the matrices $ = {ipi,...(pd) 

, and 'h = (tpi ,..., ipd)- Define the bilinear form: 

('h,$) = 7 ( 12 ) 

The Fixed Point space Sn 


{F{x)) (0) 


|f(0) ,-T<d<0 

F(a:(0),a:(-T),77) , 6> = 0 


(9) 
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Due to the S^v-symmetry of m the space where 
solutions (/)3 lie can be decomposed into the Fixed- 
point subspace where symmetry-preserving solutions 








emerge and a subspace with symmetry-breaking so¬ 
lutions, this was shown in [2. We analyze stability 
of the small-amplitude periodic solutions near Hopf 
bifurcations in the Fixed point space, these bifurca¬ 
tions satisfy assumptions ([lj)-(jnj) for it! > 1. In this 
subspace, equation ([S]) has the form: 


with general solutions: 

si(d) = cos(a;d)ci — sin(wi?)c 2 

S 2 {'&) = sin(a;i?)ci-|-cos(a;d)c 2 

ni{9) = cos(a;0)di — sin(a;0)d2 

712(d) = sin(a;0)di-I-cos(a;0)d2 


( 19 ) 


Xl = X2 

X2 = —/ia;2 + + cos 2 (j)^)xi 

-I- -I- cos2^=*=)a:iT- — + XiY sin 2^=*= 

- ^ [{xir -Xi)^ + {xir -ka;i)3cos2(/)=t] |, 

(13) 

then matrices ^0(77) and Ar{rj) in (ITII1 become: 
Mv) = Kii{-l + cos 24 >^) -n ) ’ 

^ ( Kfi{l + cos2(l)^) 0 ) ’ 

and F in ([5|) takes the form F = (/i /2)^, with fi = 
0, and /2: 

72(3:4,77) = - ^(a:ir + xi)^ sin 2(/)=‘= 

“ ^ [(^ 1 ^ “ * 1 )^ + (2^1’- + xi)^ cos2(j)^] |, 

(16) 


The coefficients ci = [cn 012]^, C2 = [c2i 022]^, di = 
[dll di2]^7 ^2 = [d2i <^22]^ can be obtained by con¬ 
sidering the boundary conditions, 

7 ^0(77)-hcos(a;T)2l^(?7) ^1^=0 

wl-I-sin(a;T)AT-(?7) J \ C2 J 

, (20) 

f ^0 (d)+cos(a;T)A^(?7) 7 ^ 1^=0 

-wl - sin(a;T)A7(?7) J V ^2 / 

the “orthonormality” condition (s, n) = /, and set¬ 
ting cii = 1 and C21 = 0, see [HKia for more details. 

It is also possible to decompose the solution Xt (d) 
to equation (0 into xt(d) = yi(t)si(i?)-k7/2(<)s2(■!?) + 
W4 (1?), where yi and 7/2 lie in the center subspace, such 
that 7 / 1 , 2(0 = (77-1,2(0),Xt( 0 )), and wt in the infinite¬ 
dimensional component subspace, thus, we have 

yi = ujy 2 + nJ{Q)F , . 

y2 = -wyi + nj (0)F 


w = A{r])wt + F{xt,r]) - n( (O)Fsi - nj{0)Fs2, 

( 22 ) 


We need the complex eigenfunctions As('d) = lajs^d), 
A*n{9) = iijjn{9), associated to the critical eigenval¬ 
ues A = iLu, and A = —lUJ with s{'9) = si{'9) -I- 152(7?) 
and n{9) = ni{9) + in2{9). These eigenfunctions 
can be computed solving the boundary value problem 
^Si,2 = =Fws2,i(7 ?), and ^ni,2 = ±0:712,1(7?), which, 
after substituting the operator ^( 77 ), becomes: 


Ao{r])si{0) ± Ar{r])si{-T) 
^ 0 ( 77 ) 52 ( 0 ) ± Arlr])s2i-T) 


and 

^0 (v)ni (0) ± ^7 {v)ni (-r) 
^0 (»7V 2 (0) ± ^7(»7)«2 (-r) 


- 0 : 52 ( 0 ) 

o:5i(0) 


(17) 


0:712(0) 

—o:ni(0) ’ 


(18) 


where 

0 ,7? g [- t , o ) 

1 d^( 7 /l(i)si( 0 ) ±7/2(<)s 2(0) ±w(t)( 0 )) ,7? = 0 . 

(23) 

The center manifold 

Following [Hdlllll], we know that w can be ap¬ 
proximated by the second-order expansion: 

W(7/I,y2)(7?) = i(hi(7?)7/J±2h2(7?)yi7/2±73(7?)y^), 

(24) 
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thus, by differentiating and substituting equa- where 
tion (| 2 ^ keeping up to second order terms, we obtain: 

w = -ujh2yl + Lo{hi - h^)yiy2 + w/i22/| + 0{y'^), ^ '= 

( 25 ) 


hi 

h2 

hs 



/ 0 -21 
C :=u}l I 0 
\ 0 21 


0 

-I 

0 


6x6 


and from equation ( 1 ^ . / \ 

dw ^ , q 

— = A{r])wt + T{w + yisi + 2/232) - (di2Si + d22S2)/2- V f2^Po J 


fi°qo \ 
f^^qo , 

fl^qo J 


( 26 ) 


From the definition of A{ri), equivalent to (TTlT) . we 
see that 


A{r])w = 


^{hiyf + 2 h 2 yiy 2 + hyl) ,de[-r,0) 

Ao(?7)w(0) + A^(?7)w(-t) ,d = 0, 


_ f di2 \ _ ( d22 \ 

■" V ^ 22^22 J ’ iv -C 22 rfl 2 ) ’ 

with general solution: 

/i(i?) = e^'^K + M cos(a;i?) + N siniuj'd). 


( 31 ) 


( 27 ) 

then, from equation (m, (Ea, (EH), and EH) , we 
can obtain the unknown coefficients h-i, /12, and /13, 
solving: 

hi = 2(-wh2 +/|°(di2Si(d) + d 22 S 2 (^ 5 ))), 

h2 = w{hi-hii)+fl'^{di2Si(d)+d22S2{'&)), 
ils = 2{u}h2+f^^{di2Sl{l}) + d22S2{'d))), 

( 28 ) where 


After substituting the general solution into ( 1 ^ we 
solve for M and N, and then from the boundary value 
problem we solving for K, 



Ph{0) + Qh{-T) = p - r, ( 33 ) 


and, 

^o(p)^i(0) + Ar{r])hi{-T) = 

2(—a;/i2(0) + /2°(di2Si(0) + ^ 2232 ( 0 ))), 

Ao{r])h2{0) + Arh2{-T) = 

w(/li( 0 ) — ^ 3 ( 0 )) + /2^(dl2Sl(0) + (i2232(0))). 


/ Ao 0 0 \ 

P := 0 Ao 0 \ -C, 

V 0 0 Ao / 

/ Ar 0 0 \ 

Q-{ 0 A, 0 , 

V 0 0 A, y 


( 34 ) 


Ao(7?)h3(0) + Ar{p)h3{-T) = 

2{u}h2{0) + /2^((il2Sl (0) + fi22S2(0))), 

( 29 ) 


and r := ( 0 /|° 0 0 ff f. 

The expressions for wi( 0 ) and wi(— t), necessary 
in ( 1231 ) . are: 


where 

2dy\ 

2dyl o' 

Equation (| 28 p is written as the inhomogeneus dif¬ 
ferential equation: 

^ = Ch + p cos(a;d) -I- q sin(a;i 9 ) ( 30 ) 


fii = ^ f 

n’ dyidy2 


and = 


wi( 0 ) = i ^(Mi -k Ki)y‘l + 2 {M 3 + 7^3)2/12/2 + + K 3 )yl 

wi(—r) = i ^(e“'"'^Ar|i -k Mi cos(wt) — Ni sm{(jjT))y‘f 

+ 2 (e“‘^'^Ar |3 -k M3 cos{u;t) - N3 sm{ujT))yiy 2 
+ (e“‘^'^Ar|5 -k Ms cos(wr) - TVs sin(a;r))2/|^ , 

( 35 ) 
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note that we only need wi(-d) since the nonlinear 
function in dUD only depends on xi; then by sub¬ 
stituting (1551) into (EU, we obtain: 

2 / 1 = W 2/2 +91(2/1,2/2;??) /ggX 

2/2= -W2/1+92(2/1,2/2;??) ’ ^ 

or 

91 = W2/2 + 0209? + 0119192 + 00292 + “3091 

+ 02i9i92 + 012919? + Oo39?, 

2/2 = —W9l + ^*2091 + ^119 i92 + ^0292 + ^309i 

+ ^219i92 + &129i92 + ^0392- 


shown in figure[TJ The black curve corresponds to sta¬ 
bility of periodic orbits near Hopf bifurcations with 
Re (A') > 0 (black curves in figured]), as we can see, 
these periodic solutions are all stable (o < 0). The 
red curve corresponds to stability of periodic orbits 
near Hopf bifurcations with Re(A') < 0 (red curves 
in figure [ 1 ]), these periodic orbits are unstable for 
9 < Ki 0.386, and stable for fj, > fj,*. Are 

also shown point A, B anc C. At points A and C, 
small amplitude periodic orbits are stable, whilst at 
point B, they are unstable. 


(37) 


In El is computed the coefficient a, which determines 
stability of the normal form (1571) . 




s¥(af + sl 


20 


-9i“ {9? + gf)-9?9f + gf9l 


20 


(38) 


i9*+'? 

where = — —--— 0 ^( 0 ,0). Periodic orbits near 
d^yid3y2 

Hopf bifurcation at the critical eigenvalue A = iw, 
will be stable if a < 0 and unstable if a > 0 . 


Numerical Results 
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Figure 1: SYMMETRY-PRESERVING BIFURCA¬ 
TION CURVES IN FIX(SAr) EOR K = 1.05. IN 
BLACK, HOPE BIEURCATIONS WITH RE(A') > 
0 AND, IN RED, HOPF BIFURCATIONS WITH 
RE(A') < 0. 


We reproduced some of the computations for the 
Hopf bifurcations in the Fixed point space for the 
case AT > 1 presented in [ 3 , because we will com¬ 
pute stability for these bifurcation curves using re¬ 
sults obtained in the previous section. In figure [T] 
(part of figure 10, in [5]) are shown the symmetry¬ 
preserving bifurcations curves in the parameter space 
( 9 ,, r) for K = 1.05, for both cases: bifurcations with 
Re (A') > 0 in black color, and with Re (A') < 0 in red 
color; we also choose three testing point for numerical 
simulation A = = (0.15,7.46), B = (0.3,11), 

and C = (0.421,7.10). 

In figure [ 2 ] is shown the coefficient a computed us¬ 
ing equation (l38l) . in the parameter space ( 9 ,t), for 
K = 1.05, related to the Hopf bifurcations curves 


In order to confirm our results, we computed 
branches of periodic solutions near the Hopf bifurca¬ 
tions points A, B, and C, using DDE-BIFTOOL [H 
[3|, along with the Floquet multipliers for a specific 
periodic solution chosen in the branch. 

A branch of periodic solutions with small am¬ 
plitude, emerging from the Hopf bifurcation point 
A = {y,,T) = (0.15,7.46) is shown in figure [Sj- (a). 
ln[3](b) it is shown the periodic solution profile psol 
at r = 7.5315. The Floquet multipliers related to 
psol are shown in figure [3]-(c). It is clear that this 
periodic solution is stable, since there is no Floquet 
multiplier outside the unity circle. 

For the point B = (/i,r) = (0.3,11), the branch of 
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0.15 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 

n 0.386 


Figure 2: COEFFICIENT a COMPUTED USING 
EQ.®, DETERMINING STABILITY OF HOPF 
BIFURGATIONS IN Fix(SAr) FOR K = 1.05, SEE 
FIGURE [TJ 

periodic solutions is shown in figure 0]- (a), in figure H]- 
(b), it is shown the profile for the periodic solution 
psol chosen at r = 11.3744, these solution is unsta¬ 
ble, because there is a Eloquet multiplier outside the 
unity circle, see figure ID- (c). 

Einally, the branch of periodic solutions near the 
Hopf bifurcation point C = (/r, r) = (0.421,7.10) is 
shown in figure [5]- (a). The periodic solution choosen 
in the branch is at r = 7.00, its profile is shown 
infigure [S]-(b). All the Eloquet multipliers shown in 
figure [S]-(c) are within the unity circle, therefore the 
solution is stable. 

Conclusions 

The reduction of the inifinite-dimensional space onto 
the center manifold in normal form, was applied to 
the Fixed point space for the Full-phase model in 
order to analyse the stability of small-amplitude pe¬ 
riodic orbits near simple Hopf bifurcations, in both 
cases, for Re (A') > 0 and Re (A') < 0, we found 
that in the first case periodic orbits which are stable 
(a < 0) can emerge, and, in the other case, unstable 
(a > 0) periodic orbits can emerge for p < p*{K), 
and stable periodic orbits for p > p*{K). The nu¬ 
merics show that the analytical results are correct. 



Figure 3: (a) BRANCH OF PERIODIC SOLU¬ 

TIONS EMERGING FROM POINT A = (/i, r) = 
(0.15,7.46197). (b) PERIODIC SOLUTION PRO- 
EILE AT = 0.15, r = 7.5315, T = 12.0364 seg, 
(POINT psol). (c) ELOQUET MULTIPLIERS EOR 
THE PERIODIC SOLUTION psol. 
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Figure 4: (a) BRANCH OF PERIODIC SOLE- Figure 5: (a) BRANCH OF PERIODIC SOLU¬ 

TIONS EMERGING EROM POINT B = (/r, r) = TIONS EMERGING EROM POINT C = (/i, r) = 
(0.3,11.001518). (b) PERIODIC SOLUTION PRO- (0.421,7.101329). (b) PERIODIC SOLUTION PRO- 
EILE AT /r = 0.3, r = 11.3744, T = 12.8506 seg, PILE AT /r = 0.421, r = 7.00, T = 8.8704 seg, 
(POINT psol). (c) ELOQUET MULTIPLIERS POR (POINT psol). (c) PLOQUET MULTIPLIERS POR 
THE PERIODIC SOLUTION psol. THE PERIODIC SOLUTION psol. 
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Although, we computed the coefficient a for a spe¬ 
cific value of K, the procedure shown is valid for all 
the parameter space where simple Hopf bifurcations 
appear. 

Finally, it is important to spotlight some points for 
further research: First, what is the nature of the so¬ 
lutions at the special point fi = at which the 

coefficient a changes sign. Second, analyze stability 
of the degenerate Hopf bifurcations at the Fixed point 
space for AT = 1, which are codimension 2, pure imag¬ 
inary eigenvalue and zero eigenvalue; and third, the 
stability of the symmetry-breaking degenerate Hopf 
bifurcations which have multiplicity N — 1. 
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